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ABSTRACT 

The increasing abundance of passive "red-sequence" galaxies since z ~ 1-2 is mirrored by a co- 
incident rise in the number of galaxies with spheroidal morphologies. In this paper, however, we 
show that in detail the correspondence between galaxy morphology and color is not perfect, providing 
insight into the physical origin of this evolution. Using the COSMOS survey, we study a signifi- 
cant population of red-sequence galaxies with disk-like morphologies. These passive disks typically 
have Sa-Sb morphological types with large bulges, but they are not confined to dense environments. 
They represent nearly one-half of all red-sequence galaxies and dominate at lower masses (< 10 10 M Q ) 
where they are increasingly disk-dominated. As a function of time, the abundance of passive disks 
with M 4 < lO n M0 increases, but not as fast as red-sequence spheroidals in the same mass range. At 
higher mass, the passive disk population has declined since z ~ 1, likely because they transform into 
spheroidals. Based on these trends, we estimate that as much as 60% of galaxies transitioning onto 
the red sequence evolve through a passive disk phase. The origin of passive disks therefore has broad 
implications for our understanding of how star formation shuts down. Because passive disks tend 
to be more bulge-dominated than their star-forming counterparts, a simple fading of blue disks does 
not fully explain their origin. We explore the strengths and weaknesses of several more sophisticated 
explanations, including environmental effects, internal stabilization, and disk regrowth during gas-rich 
mergers. While previous work has sought to explain color and morphological transformations with a 
single process, these observations open the way to new insight by highlighting the fact that galaxy 
evolution may actually proceed through several separate stages. 
Subject headings: galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

It has long been known that nearby galaxies roughly 
fall into two categories: 1) systems with prominent 
disks often exhibiting spiral structure and ongoing star 
formation, and 2) galaxies dominated by a spheroidal 
morphology — the most extreme being ellipticals — and 
harboring little if any star formation. With the advent 
of large galaxy surveys like Sloan Digital Sky Survey 
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(SDSS) it had been possible to quantify the bimodal na- 
ture of the local population both in terms of morphology 
and star formation as a function of what appears to be 
the best observationally accessible variable for separating 



the two populations, the galaxy's stellar mass (M» ; Bell 
et al.||2003MKauffmann et al.| |2003| |Baldry et aLl|2pM[ 
Balogh et al. 2004 1 . The bimodal nature of the ga laxy 



population is apparent to z ~ 1 JjBell et al 2004) and 



beyond (e.g., van Dokkum et al. 2006p , but it is clear 



that the fraction ot galaxies populating either side must 



evolve with time. Brinchmann & Ellis ( 2000 ) presented 
evidence for an increase since z ~ 1 in the global mass 
density of morphological sp heroidals at t he ex pense of 



star-forming irregulars, and Cowie et al. dT996[) showed 

the 



that increasingly if -band luminous and therefore more 
massive galaxies exhibit higher specific star formation 
rates (sSFRs) at early times. 

A number of more recent studies have confirmed and 
extended these early results, using mass and luminosity 
functions to trace the evolution since z ~ 2 of the two 
populations as partitioned b y the amount of star forma- 
tion or the rest-frame color \ Juneau et al.||2005| |Bundy 
l et al.|2006[|Borch et al.|2006||WiTlmer et al'l2T7U5[|Brown 
|et al. |2007[ |llbert et aL||2010|]M^asher et al.| | 2009|) and 
as a tuhction of morphological type (|Bundy et al.||2005| 
Pannella et al. 2006; Franccschini et al. 1120061 IScarlata 



et al. 



2007 



2010) 



2007b| lAbraham et al.||2007T 
Ibert et al.||2010||Pozzetti et a 



We 



van dor Wcl et al. 
.II2009I lOesch et al. 



The importance of environmental density has also 



2 



Bundy et al. 



been studied flBundy et al.||2006| |Cucciati et al. 2006 


Cooper et al.|12008; 


Capak et al.||2007a| van der Wei 


20081 IPannella et al.| 


2009| |Tasca et al.||2009 Bolzonella 


et al. 2009 ) . The emerging picture is one in which, begin- 



star-forming, disk-dominated systems into red, bulge- 
dominated spheroidals in which star formation has been 
quenched. As time proceeds, this transformation pro- 
cess works its way down the mass function, driving the 
apparent downward shift in the bimodality mass scale. 

Understanding the physical driver of this evolution has 
been a key challenge, leading to a search for a pro- 
cess capable of simultaneously shutting down star for- 
mation and turning disks into spheroidals. Morphologi- 
cal evolution has often been attributed to the "merger 
hypothesis" (Toomre 19771, in which galaxy mergers 
convert rot at ion- supported disks into pressure-supported 
spheroids. In a series of papers, Hopkins and collabo- 
rators demonstrate the feasibility of this idea through 
detailed simulations and link mergers with an explosive 
quasar or starburst phase that may heat or drive out cold 
gas in the spheroida l merger remnant, the reby quenching 



star formation (see Hopkins et al. |2008[ and references 
therein). The idea that mergers quench star formation, 
at least in halos above a critical mass, M cr it ~ 10 12 M Q 
(where th e cooling time is longer than the dynamical 



time, see Dekel & Birnboim 2006), has also been im- 
plemented in semi-analytic models, in which supposedly 
spheroidal merger remnants are kept "red and dead" by 
so-called "radio mode" active galactic nucleus (AGN) 



feedback ( Croton et al.|2 006; Bower et al. 2006; Cattaneo 



et al. 2006 



Somerville et al. 2008) 



Birnboim ct al. (2007 



argue against the need tor mergers and AGIN feedback, 
relying instead on the interplay between s hock heating 



and c old flows in halos near M cr j t (see also |Dekel et al 
20091. Semi-analytic models based on this framework 



have needed to include some additional quenching tied 
to bulge growth after galaxy mer gers to match the ob 
served fraction of passive galaxies ( Cattaneo et al.|2008 ) 



One way to distinguish between these scenarios and 
gain insight into the transformation mechanisms at work 
is to test the link between quenching and morphologi- 
cal evolution by studying those galaxies whose morpho- 
logical type places them on one side of the bimodal- 
ity but whose star formation identifies them with the 



Schawinski et al. 


2009; Kannappan et al. 2009; Ferreras 


et al.|2009f 


Brand et al. 


2009) and are understood to be a 



combi nation oi recent arrivals to the red sequence (Ruh 



land et al.|2009 1, bulge-dominated systems expe riencing 



triggerecTs tar formation after minor mergers (Kaviraj 
et al.|2009 ) or, at low mas ses, the formation o f new disks 



in gas-rich environments ( Kannappan et al. 2009 ) . The 
other possibility, namely quenched or passive disk galax- 
ies, is more challenging to explain because the shutdown 
of star formation must leave the disk in place. Yet ex- 
amples in the local universe do exist. Lenticular or SO 
galaxies are one example, although they typically live in 
galaxy clusters where interactions with th e intercluster 



mediu m (e.g., ram pressure stripping, see Moran et al. 
2007a) can remove the gas supply and shut down star 
formation. However, SO's have properties that do not 
depend on cluster density and a non-negligible fraction 



(~8%) of nearby SO's are isolated, suggesting that no t all 
SO's form as a re sult of living in dense environments ( van 
At slightly higher redshifts, z « 0.1, 



den Bergh 



2009) 



axis ratio distributions and Galaxy Zoo (GZ) morpholo- 
gies have demonstrated an unmistakable population of 
red, disk-dominated galaxies that appear to live in the 
field and slightly overdense environments, but are not 
found sole ly in dense regions (j B amford et al.|2009 Mas- 
ters et al |[2010| Jvan der Wei et al.|l2009p . Finally, us- 
ing the Hubble apace Telescope {TJBrj data in the Great 
Obser vatories Origins Dee p Survey (GOODS) fields at 



1, |Bundy et al. ( |2006[ ) compared color and morpho 
logically classified mass functions to argue that quench- 
ing appeared to precede spheroid formatio n, finding hints 
of a r ed disk population at z ~ 1 (see also Pozzetti et al. 
2009|. Some e xamples of passive di sks have even been 



found at z ~ 2 ( |Stockton et al.| [2008) 



The goal of this paper is to utilize the unprecedented 
size of the COSMOS field with its 1.64 deg 2 of #ST imag- 
ing to revisit the question of red, passive disks in field en- 
vironments and quantify their origin and evolution since 
z ~ 1. We argue that the behavior of this population 
presents a challenge to models that attribute morpholog- 
ical and color evolution to a single event, and may shed 
light on mechanisms that affect the galaxy population 
more broadly. This work follows on mass function studies 



logical type by Ilbert et al. 


((20T0l),|Pozzetti et al. J2009), 


Oesch et al. 


(2010), and 


Drory et ai.| (|2009 


). Ilbert et al. 



havior of morphologically early-type galaxies and those 
with passive stellar populations, but also presented evi- 
dence for an increasing fraction at high redshift of non- 
spheroidal morphologies among quenched systems — this 
is the starting point for the present study. 

Throughout this paper, we use the AB magnitude sys- 
tem and adopt a standard cosmology with Hq=70 h?o 
km s- 1 Mpc" 1 , f2 M =0.3 and Q A =0-7. 

2. OBSERVATIONS 
2.1. The COSMOS Field 

The COSMOS field (10 ,l 00 m 28.6 s , +02°12'21.0") is the 
largest contiguous area (1.64 deg 2 ) imaged by the HST 
and is almost 20 times the size of the two GOODS fields 
combined. As such, it marks an unprecedented step for- 
ward in studies of the morphologies of distant (z ~ 1) 
galaxies at high resolution. Specifically for this work, 
COSMOS provides such large samples that we can sta- 
tistically track evolving sub-samples of field galaxies de- 
fined by both rest-frame color as well as morphology as 
a function of M* and redshift from z *~ 1 to z ~ 0.1. 

A large team of scientists has assembled to collaborate 
on a wide range of topics using observations from th e 
COSMOS field (for a review, see |Scoville e!TaLp)07b ). 
The primary data sets of interest for this work are the 
isuw observations from the HST Advanced Camera for 



Surve ys (ACS; Scoville et al. 2007a Koekemoer et al 
2007 1 , which provide morphological information, and the 



ground- based optical to near- IR photometry presented in 



detail in Capak et al. ( 2007b I and updated in Ilbert et al 
(2009). As we discuss further below, the ground- based 
catalog provides the basis for high precision photometric 
redshifts and a set of rest- frame colors that separate truly 
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passive stellar populations from those that are actively 
star-forming but dusty. 
We briefly review the relevant components of the 



grou nd-based cata l og, re ferring to Capak et al. (2007b) 



and Ilbert et al. (2009) for details! The catal og in- 



cludes u'- and t'^band imaging from Mega-Prime ( Aune 



et al.|2003 Boulade et al.|2003) and was cross-referencec T 
with K 8 - band observations ( McCracken et al. I 2010| 
from the Wi de-field InfraRed Camera (WlR(JAM, [Puget 



et al. 



on 



"the 



20041), both of which are instruments 
Canada- France-Hawaii Telescope (CFHT). The major- 
ity of photometric bands in the catalog were observed 
with the Subaru T elescop e equipped with Suprime-Cam 
(iKomiyama et al.||2003t (filters Bj, Vj, g + , r+, i+, z + ). 
The catalog was selected based on detections in the Sub- 
aru i + -band and the seeing of the f i nal mo saics is better 
than l'/5 in all cases. Capak et al.| ( |2007b[ ) point-spread 
function (PSF) match all data, optimizing for photome- 
try performed in 3" diameter apertures. All of the op- 
tical bands have a 5er point-source (3" diameter aper- 
ture) depth of at least 25th magnitude (AB) with the 
-FT s -band limited to K s < 24. Given these depths, the 
final sample completeness is determined mostly by the 
ACS data because morphological classification begins to 
degrade at magnitudes fainter th an isuw ~ 24. The 



HS T ACS imaging is describ ed in Scoville et al. 
and Koekemoer et al. (2007), and is used only 



(2007a) 



phological classification in this work. We use the Zurich 
Esti mator of Structural Ty pes (ZEST) morphology cat- 
alog ( Scarlata et al.|2007a[ ) — discussed further in Section 
I3. 4J— which is based on the ACS source catalog presented 
Leauthaud et al. (20071 who estimate the ACS depth 



at isuw = 26.1 for photometry of a 1" galaxy 
2.2. The Primary Sample 

In the analysis that follows, o ur primary sample i s 
based on the --selected catalog from Capak et al. ( 2007b ) . 
We impose magnitude limits of K s < 24 and % < 2(5.5 and 
divide our sample into six redshift bins: [0.2-0.4), [0.4- 
0.6), [0.6-0.8), [0.8-1.0), [1.0-1.1), [1.1-1.2). Given the 
magnitude cuts, we estimate the mass completeness of 
each redshift bin by considering the observed magnitude 
of a maximal M„jL stellar population model with solar 
metallicity, no dust, and a r = 0.5 Gyr burst of star for- 
mation that occurred at Zf orm = 5. This exercise provides 
limits that roughly match the 80% completenes s limits 
determined when deeper samples are available (Bundy 
et al.|200"6 ), and in this case gives values of 8.8, 9.3, 9.7, 
10.0, 10.1, and 10.2 in units of log M* / M Q . 

In practice, morphological classification is not possible 
for every source down to the K s - and i-band limits. This 
means that the effective M* above which morphologi- 
cal samples are complete is higher than the full sample. 
These depend on the nature of the source, however, but 
can be easily defined with respect to the deeper K s < 24 
and i < 26.5 sample. We consider the morphological 
sample incomplete when the fraction of galaxies with a 
suc cessful ZEST morphological classification (see Section 
3.4) falls below 50% of the primary sample. We also limit 
the sample to ensure the robustness of ZEST classifica- 
tions using ellipticity distributions tests (see Section 4.2 ). 
For each of the redshift bins defined above, we adopt 
mass limits of: 9.7, 9.9, 10, 10.4, 10.5, 10.8 in units of 
log M*/M Q . After removing stars and bright AGN, the 



number of galaxies in the primary sample in each redshift 
bin is listed in Table [TJ 

2.2.1. Cosmic Variance 

While population studies using COSMOS data largely 
agree with previous work based on smaller fields (see Sec- 
tion [I]) , the sample or cosmic variance even in surveys 
covering a few square degrees represents a significant un- 
certainty t hat must be consider ed when interpreting our 
results (see Meneux et al.|[2009[ ). The largest effect is on 
the o verall nor malization of the total galaxy abundance 
(e.g., Bundy et al.||2006 l although fractional ab undance 
values are also affected. Stringer et al.| ( 2009 ) use the 
variance gleaned from multiple mock light cones based on 
the Millennium Simulation to estimate the cosmic vari- 
ance signal for a number of surveys, including COSMOS. 
From one redshift interval to another, they find typical 
abundance uncertainties of 0.2 dex. This estimate is con- 
sistent with the results in this work. 

3. METHODS 

3.1. Estimating Photometric Redshifts 

While a spectrosc opic survey, zCOS MOS, is ongoing in 
the COSMOS field flLilly et al.|[2007| ), the large samples 
required for this work necessitate the use of photometric 
redshifts. We use the catalog of photo-zs published by 



Ilbert et al. ( 2009 ) which are based on 30 bands of imag- 
ing m the COSMOS field. In addition to the observations 
mentioned in Section [2j these include UV bands from the 
Galaxy Evolution Explorer (GALEX), four mid-IR bands 
(3.6-8 /im) from the Infrared Array Camera (IRAC) on 
the Spitzer Space Telescope, J-band imaging from the 
Wide Field IR Camera (WFCAM) on the UK IR Tele- 
scope (UKIRT), and an additional 12 medium-band and 
two narrowband filters observed with Suprime-Cam on 
Subaru. This extended data set is described in detail by 
P. Capak et al. (in preparation). 
Redshift estimates were obtained with the Le Phar^\ 



code using t e mplat es drawn from Polletta et al.| (2007). 
Ilbert et al. (2009) develop a novel technique that ac- 
counts for the contribution of emission line flux (espe- 
cially important in the narrow and medium filter bands) 
to the spectral energy distributions (SEDs) of star- 
forming galaxies. This correction improves the photo- 
z precision by a factor of ~2. When comparing to 
zCOSMOS spectroscopic redshifts with 17.5 < i+ < 22.5 
Ilbert et al.| ( [2009| find a precision of a/^ z /{\ + z) ~ 
0.01 after removing an outlier fraction that accounts for 
0.6% of the comparison sample. They also compare 



to sp ec-zs obtained with Keck DEIMOS (Faber et al 
2003 1 for a sample of very faint and 24 /j,m-selected 
f 25.5 the uncertainties are estimated 
0.06. 



sources. At i 1 
at o~a.z/ (1 + z) ~ 

3.2. Estimating Stellar Masses 

We use the Bayesian code described in |Bundy et al.| 
( 2006 1 to estimate the current mass in stars of the galax- 



ies m our sample. The observed SED of each galaxy is 



compared to a g rid of 13440 models from the Bruzual 
& Chariot ( 2003[ ) (BC03) population synthesis code that 



spans a range of metallicities, star formation histories 

14 www.oamp.fr/people/arnouts/LE_PHARE.html 
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TABLE 1 
Sample Statistics 



z All Passive Star-forming Passive Early Disk Passive Late Disk Passive Extreme Disk 



0.2 < z < 0.4 


17210 


2402 


14808 


755 


380 


30 


0.4 < z < 0.6 


16324 


1744 


14580 


523 


171 


17 


0.6 < 2 < 0.8 


22253 


2821 


19432 


883 


261 


36 


0.8 < z < 1.0 


25581 


4075 


21506 


1282 


356 


81 


1.0 < z < 1.1 


6014 


480 


5534 


160 


12 


5 


1.1 < z < 1.2 


14860 


517 


14343 


165 


28 
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(parameterized as exponentials), ages, and dust content. 
No bursts are included in our models and only those mod- 
els with ages (roughly) less than the cosmic age at each 
redshift are co nsidered. We u se a Chabrier initial mass 
function (IMF; |Chabrier|2003[ ) and assume a Hubble con- 
stant of 70 km s _i Mpc \ 

At each grid point, the i^-band M*/Lpc ratios, in- 
ferred M*, and probability that the model matches the 
observed SED is stored. This probability is marginalized 
over the grid, giving the stellar mass probability distri- 
bution, the median of which is taken as the final estimate 
of M*. The width of the distribution provides the un- 
certainty which is typically 0.1-0.2 dex. This is added in 
quadrature to the -ftT s -band magnitude uncertainty to de- 
termine the final error on . Stellar mass estimates for 
galaxies with photo-zs also suffer from the uncertainty in 
luminosity distance introduced by the photo- z error and 
the pos sibility of catastrophically wrong redshift info r- 
mation ( |Bundy et al.||2005| [Kitzbichler fc White) | 2007 |. 

More broadly, any stellar mass estimate suffers po 
tential systematic erron| 15 
in stellar population mod"e 



by recent work showing how this approach separates 
dusty star-formers from tru l y passive populations (e.g., 
Pozzetti fc Mannuccil [20001 ILabbe et al.| [20051 |Wuyts| 



po- 

from uncertainties inherent 
ing and various required as- 



sumptions, such as the form of the IMF. Several papers 
have stressed the importance of treating thermally pul- 
sating asymptotic gia nt branch stars (TP-AGBs, e.g., 
Maraston et al.|[2006[) an element that is missing in the 
Bruzual & Chariot ( 2003]) models. The recent and thor- 
ough investigation o f population synthesis modeling by 
Conroy et al. ( 2009), however, argues th at M/L ratios 
estimated from Bruzual & Chariot ( 2003 ) are largely re- 
sistant to the uncertain contribution from TP-AGBs as 
well as other limitations of curre nt models. A further test 
is provided in |Drory et al.| ( |2009[) who compare mass func- 
tions computed using BC03 and BC07 stellar population 
codes and conclude that ^20% differences in abundances 
can be introduced, especially for star-forming galaxies. It 
is also important to recognize that estimates may be 
affected by unrecognized systematic uncertainties at the 
0.1-0.2 dex level. 

3.3. Excluding Dusty Red-sequence Galaxies 

Because our goal in this paper is to study disk galaxies 
that have stopped forming stars, in this section we dis- 
cuss the rest-frame galaxy colors we use to identify pas- 
sive stellar populations. Previous work has often relied 
on a single rest-frame UV-optical color to select passive 
galaxies, but this technique leads t o possible contam ina- 
tion from dusty star- formers (e.g., 
this work, we make use of color-co 



Yan et al.||2006[ ). In 
or diagrams that in- 
clude rest-frame UV, optical, and near-IR colors, inspired 

15 These are not included in our final error estimates. 



et al. [20071 [Williams et al.||2009p 

The utility of the rest-frame ITV-optical-near-IR "bi- 
color" diagram for the COSMOS data set is illustrated 
in Figure [l] where we show interpolated rest-frame colors 
in various redshift bins using the following filter bands: 
the GALEX near-UV band (NUV), the Suprime-Cam r+ 
band which we label as R, and a generic J-band filteip'l 
Rest-frame N UV- and .R-band ab solute magnitudes were 
estimated by Ilbert et al. ( 2009 1 during the SED fitting- 
necessary for measuring photo-zs. To derive the rest- 
frame J-band absolute magnitudes needed for the NUV- 
R-J diagram, we used t he Kcorrect software package 
( Blanton fc Roweis]|2007|_ and the full suite of photome- 
try described in Section [51 The broad wavelength cover- 
age of COSMOS ensures that the flux corresponding to a 
given rest-frame filter is well sample d. A similar diagram 
was discussed in Ilbert et al.| ( |2010 ) . 

The upper- left portion of each redshift panel in Figure 
[I] exhibits a relatively tight clump of galaxies that we 
associate with passive stellar populations. A sequence of 
star-forming galaxies is evident below and to the right 
of the clump in each panel. As a test, we have colored 
in red those galaxies with SEDs best fit by early "spec- 
tral" type templates, as determined by the photo- z code 
(this cate gory includes types 1-8 in the COSMOS photo- 
z fits, see Ilbert et al. 2009). Highly star-forming tem- 
plates are shown in blue. At most redshifts, the early- 
type SED designation overlaps well with the associated 
passive clump. Note that the location and orientation 
of the passive clump and star-forming sequence evolves 
somewhat given our choice of filters. In what follows, we 
will isolate passive clump galaxies (also referred to here- 
after as the "red sequence") using a joint color cut that 
evolves accordingly and is indicated by the two lines in 
Figure [T] Recognizing that such color cuts are somewhat 
arbitrary, our goal is to evenly split the now meandering 
"green valley" that separates the two populations. We 
specify our thresholds in terms of the given filters be- 
cause they correspond well to our observations and do a 
good job of making the passive clump distinct, although 
we acknowledge that a more general treatment would be 
valuable in future work. Still, we emphasize that adjust- 
ing the details of the cuts defined below has little effect 
on the analysis that follows. 

16 As this work was completed, J-band observations of the COS- 
MOS field were obtained using the WFCAM on UKIRT which has 
a different filter response function than the generic J-band filter 
used here. 
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The horizontal NUV— R color cut is tuned to the usual 
lumin osity dependent bimodality 



(e. 



Willmer et al. 



2006) in this color and is defined as follows 



NUV — R > 4.2(1 + z) 



0.2(M x + 20), (1) 



where Mr- is the rest-frame absolute -fT s -band magni- 
tude. Passive galaxies must also satisfy the following 
diagonal cut in Figure [TJ 



NUV - R > CiO) + C 2 (z)(R - J), 



(2) 



where the constants, C\ and C2, have been chosen by 
inspection in redshift bins of width 0.2 with central val- 
ues of z = [0.30,0.50,0.70,0.85,0.95,1.10]. They are 
given by d{z) = [4.4,4.2,4.0,3.9,3.8,4.2] and C 2 {z) = 
[2.41,2.41,2.5,2.6,3.0,3.7]. At z < 0.6 the combination 
of the color cuts above select 10%-20% more passive 
galaxies than those with early-type SED templates, many 
of them scattered between star-forming and passive se- 
quences. At higher redshifts, however, the SED method 
appears to miss genuine galaxies in the passive clump, 
as many as 20% at z ~ 1. We note that our conclusions 
regarding passive disk galaxies are unchanged if we use 
the SED fits to define the passive population. 

Finally, we have checked for obscured star formation 
among the passive population as indicated by MIPS 
24/im emission. We note that a significant fraction of 
24/im emission, especially am ong the passive p opulation, 
might arise from AGN (e.g., Fiore et al.||2008[) . In ad- 
dition, among star-forming galaxies, the implied SFR 
from 24/im emission must be added to that from the 
UV flux which will dominate the ou tput luminosity from 
star formation in our sample (e.g., Noeske et al. 2007). 
As a result, our constraints on hidden star formation 
in passive galaxies compared to star-formers are upper 
limits. Using the morphological classifications defined 
below (grouping "early" and "late" disks together), our 
first test compares 3c MIPS detection rates of passive 
versus star-forming disk galaxies. For z > 1, no passive 
disks are detected with MIPS while the detection rate 
among star-forming disks is roughly 10%. In our low- 
est redshift bin, 0.2 < z < 0.4, 1% of passive disks are 
detected compared to 10% of star-formers. At redshifts 
in between, the highest detection rates of passive disks 
occur at 0.4 < z < 0.6, where 5% are detected compared 
to 8% for star-forming disks. Considering the full red- 
shift range, these numbers confirm that the color-color 
methodology described above does not lead to signifi- 
cant contamination from hidden star formation. How- 
ever, given the small numbers of MIPS detections, it is 
not possible to compare 24/im flux as a function of M* 
or put constraints on the relative amount of SFR hidden 
among passive population but below the MIPS detection 
limits. 

For this, we perform a stacking analysis of the MIPS 
data in bins of stellar mass (of width 0.5 dex) and in the 
redshift intervals used throughout (the min imum num- 
ber o f galaxies per bin is roughly 20). Using ( |Rieke et al. 
|2009[ ), we estimate mean sSFR given the stacked 24/xm 
emission. Across our ~30 stacked bins, the average pas- 
sive sSFR is roughly an order of magnitude less than 
the star-forming sSFR. In two bins this difference ap- 
proaches a factor of 2, and in several other bins, there is 



no stacked MIPS detection for passive galaxies. We note 
that the IR/SFR calibrations become uncertain at low 
24/im luminosities below log L24/ Lq < 8.5 and several 
of our bins reach this limit. With this caveat and again 
ignoring contributions from AGN, our stacking analysis 
provides a rough upper limit to the possible amount of 
hidden star formation (SF) in the average passive galaxy 
that corresponds to growth rates of 3%-5% in M* per 
Gyr. 

3.4. Morphological Classification 

Understanding the morphological makeup of the 
galaxy population is crucial to our goals in this paper. 
Given the large samples (typically more than 10 5 galax- 
ies) made possible by the COSMOS data, automated 
morphological classifiers figure prominently. It is well 
known, how ever, that techniq ues such as CAS 



differ from each other. Indeed, Ilbert et al 



rphqlog 



Con- 



known, how ever, that techniques such as tAb (e.g. , Con- 
selice|2003| or Gini/M20 (see [Abraham et al. \2007\ often 
disagree when compared to visual morphologies and can 



2010J) use two 
described in 

Abraham et al.| (|2007|) and a CAS-like technique used 



different classifiers — th e Gini/M20 metho 



in jCassata et al. ( 2007 ) — to gauge the systematic uncer- 
tainty of broad categories such as "spheroidal" or "disk- 
like." Their analysis of the evolving distribution of mor- 
phologies in COSMOS from z ~ 1 serves as the basis for 
this work. We should also not e that "wiki-morpho logies" 
of the kind pioneered by GZ ( |Lintott et aT|2 008) are an 
alternative way to obtain visual morphologies tor large 
samples. 

In this work, we are specifically interested in disk com- 
ponents living in passive galaxies. As such, the exact 
definitions of morphological type are less important than 
confidence in the ability of the chosen classifier to find 
disk-like structure even if embedded in bulge-dominated 
galaxies and, ideally, to parametrize the prominence of 
disks compared to an accompanying bulge component. 
The ZEST automated classifier satisfies these require- 
ments and has also been trained to match visual classifi- 
cations of local surveys. Our strategy is to adopt the 
ZEST classification in the COSMOS sample and test 
it both through visual inspection of COSMOS galax- 
ies and an independent check of our results using visual 
morphologies determ ined in GOODS-N and published in 



Bundy et al. (2005) 



3.4.1. The ZEST Automated Classifier 

ZEST combines the power of a principle component 
(PC) analysis of nonparametric measures of galaxy struc- 
ture with inf ormation from a param etric fit. Details are 



presented in Scarlata et al. (2007a) but we review the 
salient features here. The light concentration, asymme- 
try, Gini coefficient, second-order moment of the bright- 
est 20% of the flux, M20, and the ellipticity are input 
to a PC analysis which delivers eigenvectors (linear com- 
binations of the input parameters) that maximize the 
observed sample variation, thereby allowing for a com- 
pression of the original parameter space. Only the first 
three PCs, PCi, PC2, and PC3, are retained for further 
analysis because they contain the bulk of the sample vari- 
ance. The PC1-PC2-PC3 space defined by a training 
sample of 56,000 isi4W < 24 galaxies was divided into 
a three-dimensional grid. Every galaxy living in a given 




Fig. 1. — Rest-frame (NUV — R) vs. (R — J) color— color diagrams in three redshift bins selected to represent the full z-range of our sample 
(0.2 < z < 1.2). As shown by |Williams et al7| ( |2009"] | , the addition of a near-IR, filter helps break the degeneracy between dusty star-forming 
galaxies (red in NUV— R and red in K — J ) and truly passive populations, which form a tight locus that is red in NUV— R but blue in 
R—J. The lines indicate the adopted color cuts. Red and blue colors indicate early and highly star-forming SED types, respectively, as 
identified by the photo- 2 fitting code of |Ilbert et al,| ( |2009} . 

unit cube was visually inspected so that a broad mor- 
phological type could be assigned to that cube. A value 
of Tzest = 1 corresponds to "elliptical," Tzest = 2 
to "disk," and Tzest = 3 to "irregular." In the case 
of Izest = 2 (disk) cubes, the average value (n) from 
single-Sersi c fits fo r galaxies with isi4w < 22.5 from |Sar-"1 
gent et al. (2007) was used to establish a "bulge" pa- 
rameter, -Bzest, for each cube which was then applied 
to fainter magnitudes. Higher values of -Bzest indicate 
smaller bulges and larger disks, with -Bzest = [0,1, 2, 3] 
corresponding to Sersic n ranges of (n) > 2.5, 1.25 < 
(n) < 2.5,0.75 < (n) < 1.25,0 < (n) < 0.75. In the 
analysis that follows, we consider three ZEST morpho- 
logical types: 



• "Spheroidal" includes Tzest = 1 ellipticals as well 
as highly bulge-dominated SO-like systems with 
Tzest = 2 and B ZEST = «») > 2.5). 

• "Early disks" are likely to host significant bulges 
and have T Z est = 2 and -Bzest = 1 (1.25 < (n) < 
2.5). 

• "Late disks" likely have smaller bulges with 
Tzest = 2 and B ZEST = 2 (0.75 < (n) < 1.25). 

• "Extreme disks" have little to no bulge component 
with Tzest = 2 and B Z est = 3 (0 < (n) < 0.75). 

The ZEST morphologies were previously tested using a 



1996). 



de Vaucouleurs et al. 111991; Frci et al. 


in Figure 6 ot Scarlata et aI7 


(2007a 


), 



morphologies. Of relevance to the current work, galax- 
ies with Tzest = 2 (disks) and -Bzest > are always 
associated with types later than SO with no contamina- 
tion from ellipticals. Of disks with Bzest = 1, 43% have 
published types between SO and Sab and 57% have types 
Sb-Scd. Of disks with -Bzest = 2, 5% are classified as 
SO-Sab, 45% as Sb-Scd, and 50% as Sd type and later. 
At faint magnitudes the ZEST morpholo gies begin 



( 2007a ) . This can be quantified by degrading the signal- 
to-noise ratio (S /N) of bright galaxies for which the mor- 
phology is well determined. Down to iaiAW — 22.5 
fewer than 10% experience any change in either Tzest 
or -Bzest and those that do mostly lie near classification 
boundaries. Even when degraded to isuw — 24, the 
effective limit for morphological classification, the frac- 
tion identified with their original morphological type is 
65%. Because we will be interested in trends between 
morphology, mass, and redshift, it is important to note 
that at faint magnit udes, the ZEST morph ologies are bi- 
ased to later types. Scarlata et al. (2007a) estimate that 
at most 30% of early-type galaxies could be misclassi- 
fied as disk or irregular type at the faintest magnitudes. 
We will show that this potential bias is smaller than the 
detected amount by which disks outnumber spheroidals 
at the faint ends of our sample. Furthermore, imposing 
brighter magnitude limits (e.g., isnw < 22.5) to reduce 
the potential for poor classifications does not alter our 
basic conclusions. 

3.5. Estimating Comoving Number Densities 

We compute mass functions in this work using the 
Vraex technique (Schmidt 1968). We weight galaxies by 
the maximum volume in which they would be detected 
within the -ftT s -band and i + -band limits in a given redshift 
interval. In practice, our magnitude limits (K < 24 and 
i + < 26.5) impose restrictions only in the highest red- 
shift bins (z > 1) and at masses below the limit where 
morphological classification becomes difficult. For each 
host galaxy i in the redshift interval j, the value of V^ ax 
is given by the minimum redshift at which the galaxy 
would drop out of the sample 



V 1 



Zhi E h 



dQ——dz, 
dz 



(3) 



where d£l is the solid angle subtended by the survey area, 
and dV/dz is the comoving volume element. The redshift 
limits are given as, 



to break down, as discussed further in Scarlata et al 



^high 
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(4) 
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(5) 



where the redshift interval, j, is defined by [z 3 min , zj,ax] 
and z~L and zf refer to the redshift at which the 

« lim * lim 

galaxy would still be detected below the if s -band and 
i + -band limits. We use the best-fit SED template as de- 
termined by the stellar mass estimator to calculate these 
values, thereby accounting for the ^-corrections neces- 
sary to compute V max values (no evolutionary correction 
is applied). 



4.1. 



4. RESULTS 

Morphological Makeup of the Red Sequence 

Our main results are shown in Figures [2] and [3] which 
present the absolute mass function and fractionalcontri- 
bution of passive galaxies (which we refer to as the "red 
sequence") divided by ZEST morphological type into 
spheroidals, early disks, and late disks. The mass distri- 
butions of all red-sequence galaxies (red dashed line) , all 
star- forming galaxies (blue dash-dotted line) , and the full 
sample (black solid line) are also indicated. The statisti- 
cal uncertainties are encoded by the width of the shaded 
regions. Th e sam ple (or cosmic) variance in COSMOS 
(see Section 2.2.1[ ) is estimated at 0.2 dex and is likely 



responsible tor the high apparent abundance values (es- 



feature seen in previous work ( 


Ilbert et al. 20101 iDrory 


et al. 


2009| |Pannella et al. 20( 


W Fozzetti et al. (2009} . 



do not plot them and for clarity we omit the morpho- 
logical mass functions of star-forming galaxies. Though 
we do not show it, it is important to note that at low 
masses star-forming disks of all three types far outnum- 
ber passive disks, while at the highest masses they are 
roughly equal in absolute numbers. The 50% complete- 
ness of the ZEST morphological sample is indicated with 
vertical dotted lines in Figure [2] while data points below 
this limit are simply omitted in Figure [3| These lim- 
its are nearly identical to the 50% completeness limits 
that would be obtained imposing a brightness limit of 
*814VK < 24. 

Figures [2] and [3] show that galaxies on the red sequence, 
even after removing those contaminated by dust, are re- 
markably diverse in their morphologies. This may have 
been missed in previous data sets that were limited to 
the highest masses where spheroidals do indeed domi- 
nate the red sequence. Thanks to the greater depth of the 
COSMOS photometry, the ZEST morphologies reveal a 
significant population of red sequence galaxies at lower 
masses that appear to harbor disk components. These 
become prominent below a stellar mass scale that evolves 
downward with time from l-2xlO n M0 at z « 1.2 to 
«1O 1O M at z w 0.3. Considering the full red sequence 
(defined by the mass range over which the fraction of 
passive galaxies is / rc d > 10%) we find that at every red- 
shift, passive disks (early + late) are identified in as much 
as half of red-sequence galaxies and become dominant at 
the low-mass end. This increase in passive disks is robust 
to the maximum combined biases that could result from 
morphological classifications as well as cosmic variance. 
Dividing disks further into "early" and "late" depend- 
ing on whether the average Sersic (n) values are greater 
than or less than 1.25, we see that disk galaxies with 



prominent disk components (lower (n)) are more com- 
mon at lower masses (<10 10 M@) becoming comparable 
to the more bulge-dominated early disks with higher (n) 
values. On the other hand, these early disks are more fre- 
quent at higher masses and there is a redshift-dependent 
mass range (near ~10 n -MQ) where they are roughly as 
abundant as passive spheroidals. This general behavior 
seems to mirror the broader pattern seen when the to- 
tal population, regardles s of rest-frame color is split by 
morp hological type (e.g., Bundy et al.||2005 Ilbert et al. 
2009| . 

The evolution of morphologies on the red sequence is 
most easily seen in Figure [3] where the fraction of passive 
galaxies with a given morphology is plotted as a function 
of redshift in four bins. While spheroidals dominate 
the highest masses of the red sequence at all times, at 
lower masses they start off as a small fraction of passive 
galaxies at early epochs and rise in importance with time. 
Meanwhile, the fraction of early and late passive disks 
tends to decline with time. As we discuss below, it seems 
plausible from Figure [3] that the passive disks' fraction 
declines because they transform into passive spheroidals, 
but Figure [2] shows that their absolute numbers evolve 
less strongly. It is therefore likely that some star-forming 
systems (including star-forming spheroidals) transform 
directly into passive spheroidals. This is supported by 
Figure [2] which shows that at every redshift, there is an 
ample supply of star-forming galaxies which may later 
transform into passive spheroidals. 

4.2. Verification of the Passive Disk Population 

Figures [2] and [3] reveal a surprisingly abundant pop- 
ulation of passive galaxies that appear to contain disk 
components. As we argue in Section [5j these galaxies 
offer valuable clues to the mechanisms driving galaxy 
evolution. We have visually inspected early and late pas- 
sive disks identified by ZEST in COSMOS and present 
random examples drawn from specific redshift and M„ 
ranges in Figures [3] and [5j Early disks tend to be bulge- 
dominated as expected, with many examples that might 
be classified as Sa/SO galaxies. These galaxies often 
show rounded central regions surrounded by lower sur- 
face brightness material that we attribute to the disk 
component responsible for the relatively low Sersic values 
(1.25 < (n) < 2.5). Note that our choice of ZEST classi- 
fications groups objects exhibiting even higher n values 
((n) > 2.5) in the spheroidals category. The late disks 
(Figure [5]) tend to be more elongated with smaller bulges 
and comparably brighter disk components. In either cat- 
egory and across the full redshift range, there is little 
evidence for the presence of significant amounts of dust 
and only a small degree of spiral structure. Compared to 
the clumpier appearance of star-forming disks with spiral 
arms and associated star forming regions, the morphol- 
ogy of passive disk galaxies is consistent with the notion 
that star formation has been largely shut down in these 
systems. 

While direct inspection supports the automated ZEST 
clas sifica tions, systematic biases could remain. In Sec- 
tion [3T4J we discuss how at faint magnitudes ZEST may 
misclassify up to 30% of spheroidals as disks. This po- 
tential bias is illustrated by the systematic error bar in 
Figure [2] which shows that passive disks would still dom- 
inate the low-mass end, even if the maximum 30% cor- 
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Fig. 2. — In six redshift intervals, mass functions of galaxies on the red sequence divided by morphological type. The total mass function 
is shown by the solid black line (the lowest redshift mass function is repeated as a gray line in all panels) while the abundance of passive 
galaxies (of all types) is shown as the dashed red line and that of star-formers with the blue dash-dotted line. The passive (red-sequence) 
mass functions are further divided by morphological type as shown by the shaded regions which indicate the la uncertainties. Passive 
E/SOs have red shading with asterisks, passive early disks (1.25 < (n) < 2.5) have purple shading with triangle symbols, and passive late 
disks (0.75 < (n) < 1.25) have orange shading with X symbols. The isolated error bars indicate the systematic uncertainty associated 
with a 30% upward correction to the number of ZEST spheroidals potentially misclassified at the faintest magnitudes. The contributions 
from irregulars and passive galaxies with extreme disks (which are very rare) are not shown for clarity. Cosmic variance tends to affect the 
abundances of sub-populations to a similar degree and introduces differences between redshift panels that can be as large as a factor of 4 
but are typically closer to 60%. 



rection were applied. We have explored this bias fur- 
ther by comparing the ellipticity distributions of passive 
and star-forming disks. At the mass completeness lim- 
its of our analysis, the ellipticities of passive and star- 
forming early disks are similar, but late passive disks are 
weighted toward lower ellipticities. This occurs at a level 
that is consistent with the expected 30% contamination 
of the late disk sample by misclassified early disks and 
spheroidals. Thus, at the mass completeness limits and 
below, the passive late disk abundance should be viewed 
with caution given that as much as 30% could be mis- 
classified. 

We emphasize, however, that the basic trends in our 
results are robust because they are still present when we 
limit our sample to bright magnitudes (isuw < 22.5) 
where the classifications are hi ghly secure. In abs olute 
terms, we see by comparing to Ilbert et al.| ( |2010[) that 



disks as spheroidals. This is expected given that broad 
divisions are somewhat arbitrary and that our goal is 
to identify evidence for disk components even in bulge- 
dominated galaxies that might otherwise be grouped 
with spheroidals. The galaxies with disk morphologies 
in our sample are associated with average Sersic indices 
of (n) < 2.5 for the early disks and (n) < 1.25 for late 
disks. A value of n — 2.5 h as often been used t o divide 
disks from spheroidals (e.g, Barden et al. 2005). How- 



other automated methods classiiy some fraction ot ZEST 



ever, we should note that surface brightness dimming 
at high redshifts may make some disk galaxies appear 
more bulge-dominated, although this works in the oppo- 
site sense as the estima ted bias at faint magnitudes from 
Scarlata et al. (2007b). So, until verification through 
bulge-to-disk decompositions is performed, a conserva- 
tive interpretation would be to view the abundance of 
ZEST disks as an upper limit. 
A more specific issue is a possible bias from inclination. 
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Fig. 3. — In 4 M, intervals, the fractional contribution to all 
red sequence (passive) galaxies from E/SO's, early-type disks, and 
late- type disks as a function of redshift. The contribution from 
irregulars is not shown for clarity. Only data points above the 50% 
completeness limit for obtaining successful morphological classifi- 
cation are plotted. 

We expect any type of disk component to be more eas- 
ily identified in edge-on systems where natter isophotes 
are easier to distinguish. Of particular concern is the 
possibility that dust contamination in edge-on systems 
could hide residual star formation given that the inte- 
grated colors could be dominated by a quiescent bulge. 
To check this, we have recomputed the mass functions 
of passive galaxies with the requirement that all galaxies 
have a ZEST ellipticity class of 0, corresponding to com- 
pletely face-on systems. This reduces the sample size 
by ^80%. While the mass functions shift a bit under 
this assumption, the essential behavior captured in Fig- 
ure[3]is unchanged, with the same relative abundance and 
mass dependence of the various morphological types pre- 
served. We can even compute mass functions requiring 
significant elongation (ellipticity class 3). While this dis- 
favors spheroidals because they are intrinsically rounder 
and boosts late disks which are intrinsically flatter, the 
basic mass and redshift trends are still unchanged. We 
can therefore rule out an inclination bias in our ability 
to identify passive disks. 

Another concern is the morphological ^-correction 
which could mimic the rising abundance of passive disks 
at early times if, as expected, galaxies appear more disk- 
dominated at bluer rest-frame wavelengths. Over the 
redshift range of the sample (0.2 < z < 1.2) the ACS 



*8i4W filter corresponds to a rest-frame wavelength range 
of 4000-6700 A. It has been shown that morphological K- 
corre ctions within this range are small fo r most galaxies 
(e.g., |Lotz et al.|2004]|Cassata et al.|2005| , but with only 
one b and of ACS data, it is hard to test this effect di- 



rectly. Capak et al. (2007a) use a subsample of COSMOS 
galaxies to show that automated estimators are generally 
robust to ^-corrections, but to address this problem fur- 
ther and check our ZEST-based results more generally, 
we use the visua lly classified morpholog ies in the north- 
ern GOODS- N (Giavalisco et al.||2004[ ) as presented in 



Bundy et al. 



(2005). These observations have a compa- 
rable depth as the COSMOS isuw imaging, but include 
4 HST/ ACS bands (BViz), enabling morphological clas- 
sifications robust to potential if- correct ions. 

We identify quiescent stellar populations using simi- 
lar color-color cuts for the GOODS data as for COS- 
MOS and do indeed find a population of passive disks 
with visual morphological types of Sa-Scd. Figure [6] 
shows several examples. While the small size of the 
GOODS comparison sample prevents a detailed compari- 
son, it nonetheless shows similar evidence for an increas- 
ing fraction of passive disks as a function of redshift. 
For 0.4 < z < 0.7 and above the completeness limit of 
log M*/Mq w 10.2 we find <20% of the quiescent pop- 
ulation is identified as disk-like. At 1.0 < z < 1.1 and 
logM*/M Q > 10.7, the fraction rises to roughly 35%. 
While the statistical uncertainty is large, the sense of 
this trend reinforces the results of the larger COSMOS 
sample. It should be noted that biases among visual mor- 
phologies favor early type classifications at faint magni- 
tudes and high redshifts due to sur face brightness dim- 
ming (e.g., |Brinchmann et al. 19981. We also note that 
the GOODS objects show no apparent bias with inclina- 
tion, in agreement with the inclination tests above. The 
use of color information in refining the GOODS mor- 
phologies indicates that the evolutionary signal is not 
likely to result from If -corrections and the fact that pas- 
sive disks identified visually exhibit similar abundances 
and evolution lends support to the ZEST classifications. 

5. DISCUSSION 

We have used automated Jsi4iy-band ACS morpholo- 
gies in the COSMOS survey supported by visual morpho- 
logical classifications of BViz ACS imaging in GOODS- 
N to argue for the existence of a significant population 
of passive, bulge-dominated galaxies hosting disk com- 
ponents. These galaxies were more abundant in the past 
(z ~ 1) and are seen to be more common among galax- 
ies with lower at all redshifts probed. In fact, at 
the low-mass end they represent as much as 70%-90% 
of the red sequence, even after dusty star-forming red 
galaxies are removed. At these masses, passive disks 
account for about 10% of all galaxies. The fact that 
their abundance depends on both M* and redshift, and 
that these dependences are influenced by the strength of 
the disk component argues that passive disks hold clues 
to understanding both the quenching of star formation 
and morphological evolution within the galaxy popula- 
tion more broadly. We begin our discussion with an as- 
sessment of how important passive disks are to galaxy 
evolution, finding that as many as 60% of galaxies tran- 
sitioning onto the red sequence may evolve through a 
passive disk phase. We then discuss possible formation 
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Fig. 4. — Examples of passive early disks identified using ZEST in the COSMOS data. The left-hand panel shows galaxies with 
0.2 < z < 0.4 and 9.7 < log M*/M < 10.2 while the right-hand panel corresponds to 0.6 < z < 0.8 and 10.2 < IogM*/M Q < 10.7. This 
class features bulge-dominated systems with average single-Sersic n values of 1.25 < (n) < 2.5. The galaxies displayed were chosen at 
random from the subsample. Their redshifts are indicated in the upper left of each postage stamp, the M t values in logM0 units in the 
lower left. Each panel is 3" across, corresponding to roughly 12 kpc at z ss 0.3 and 21 kpc at z ss 0.7. 
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Fig. 5. — Examples of passive late disks with the two panels corresponding to the same mass and redshift range as in Figure [4] These 
objects have average Sersic n values of 0.75 < (n) < 1.25 and are more disk-dominated than the early disks. The postage stamp sizes and 
labeling is the same as in Figure [4] 
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Fig. 6. — Viz color images of passive galaxies in GOODS-N with visually classified disk morphologies. The upper strip of eight postage 
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indicated in the upper left, logM„/Mo in the lower left. 
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mechanisms, first ruling out simple fading of star-forming 
disks and strong environmental processes and then mov- 
ing on to discuss disks regrown in mergers and the sug- 
gestion that quenching is initiated by internal structural 
changes. While no current model fully reproduces the 
observations at present, our constraints will help refine 
proposed quenching mechanisms and indicate that mor- 
phological and color transformations in galaxies proceed 
in multiple and sometimes separate stages. 

5.1. The Importance of Passive Disks as a Phase of 
Evolution 

How significant are passive disks in the evolutionary 
history of galaxies? Ideally, we could answer this ques- 
tion by tracking the mass-dependent flow of galaxies 
transitioning through different evoluti onary phases, but 
unfortunately, as discussed in Section 2.2.1 the strong 
cosmic variance uncertainties in the single COSMOS field 
prevent a detailed treatment. Broadly speaking, Fig- 
ures [2] and [3] in conjunction with previous work suggest 
that passive disks may represent a significant and rel- 
atively common phase of evolution. While the red se- 



quence abundance grow s with time (e.g., Bundy et al 



2006[ |Borch et al 
at lower masses w 



2006 ), passive disks always dominate 
rich is also where the newest arrivals 



to the red sequence can be found (e.g., |Treu et al.|2 006) 
This suggests that a large fraction or these quenching 
galaxies experience an initial passive disk phase. This is 
supported by Figure [2] which shows that, as with red- 
sequence spheroidals, the absolute number density of 
passive disks below ~10 n MQ increases with time, even 
though, compared to all passive galaxies in this mass 
range, their fraction decreases with time. Choosing a 
mass between 10 10 and 10 11 Af Q and comparing the high- 
est and lowest redshifts sampled, we see that ^20%-40% 
of new red sequence galaxies formed over this interval 
can be accounted for by the increase in the number of 
passive disks. At masses above ~10 11 M Q , evolution in 
the absolute abundance of passive disks is more difficult 
to discern but, as with their fraction, likely declines by 
20%-50% (compare the two lower redshift bins in Fig- 
ure [2] and note the lack of high-mass examples at late 
times). This reinforces the visual impression of Figure 

f namely that passive disks evolve into spheroidals. In- 
5ed, as we discuss below, one expects morphological 
evolution of passive disks to be enhanced compared to 
their star-forming counterparts because lower gas frac- 
tions (/gas) enable more minor an d, hence, more frequen t 
mergers to drive transformations ( |Hopkins et al.|2~009a [). 
If we make the assumption that all passiv e disks trans- 
form into sp heroidals within 1-2 Gyr (e.g., van Dokkum 



et al. 2005), the fraction of new arrivals to the red se- 



quence transitioning through a disk phase must increase 
to compensate, and could be as high as 609?[[^J Pas- 
sive disk galaxies may therefore represent a fairly com- 
mon phenomenon and highlight the importance of fur- 

17 For a mass between 10 10 and 10 11 Mq this percentage can be 
estimated by assuming that all of the passive disks at z Ri 1 become 
passive spheroidals by z ss 0.7 (an interval of roughly 1.5 Gyr). The 
remaining spheroidals that have formed between these two redshifts 
therefore originate from star-forming galaxies. This accounts for 
~40% of the total growth of the red sequence abundance over this 
interval. Thus, ~60% of red sequence growth must take the form 
of passive disks. 
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Fig. 7. — Concentration of passive vs. blue disks of compa- 
rable morphologies in the 0.2 < z < 0.4 redshift range with 
9.7 < logM/AfQ < 10.2. Passive disks in the same morpholog- 
ical class show evidence of higher concentrations. The ACS igi4w 
band samples the restframe at ~6200A, thus avoiding potential 
classification biases owing to star-forming regions. 

ther morphological evolution along the red sequence. 

5.2. Can Passive Disks Form From Faded Star-forming 

Disks? 

At all times, there are as many if not more star-forming 
disks than passive disks, so it is easy to imagine that 
passive disks may be drawn from the star-forming pop- 
ulation with one important caveat: there are almost no 
passive galaxies with "extreme disk" (Sersic (n) < 0.75) 
morphologies. Such extreme disks are common in large 
numbers among star-forming galaxies at all masses and 
even dominate at the lowest masses where passive disks 
are most abundant. They would be easily detected if 
they existed in the passive population. Therefore, the 
transformation to the red sequenc e either leads to the for- 
mation a bulg e component (e.g., Dommguez-Palmero & 
Balcells|2 009)' or requires the existence of one. In the lat- 
ter case, star-forming and passive disks with similar bulge 
components should differ only in color but not in their 
intrinsic structure. We will postpone a detailed com- 
parison of bulge/disk decompositions, scale lengths, and 
M/L evolution to future work and focus here on a sim- 
ple comparison of the concentration (C = Rso%/Iho%) 
of passive and star-forming early disks (Figure Yf\ . We 
consider the nearest redshift bin 0.2 < z < 0.4 so that 
the ACS i8i4M/-band probes the reddest possible wave- 
lengths (least affected by star formation), and consider 
the mass range of 9.7 < log M*/M < 10.2. Passive 
disks appear to be moderately more concentrated with a 
peak at C « 3.9 compared to blue disks with a peak of 
C ~ 3.6 (a Kolmogorov-Smirnov probability test shows 
that the probability of 0.02 for the distributions being 
drawn from the same sample). Even after dividing the 
star-forming disks into low and high SFR subsamples 
(based on their NUV luminosities) and comparing just 
the low SFR star-formers to passive disks, this difference 
in concentration remains. It also holds to z 0.9 despite 



12 



Bundy et al. 



rest-frame morphologies that sample bluer wavelengths. 
Still, this signature could largely be driven by differences 
in M/L gradients between star-forming and passive sys- 
tems. Given this tentative result and the stronger re- 
quirement that passive disks be more bulge-dominated 
than star-forming disks, it seems unlikely that passive 
disks form simply via the quenching and fading of star 
formation disk galaxies. 

5.3. The Role of Environment 

Could passive disks be a product of dense environ- 
ments? Their morphologies are sometimes reminiscent 
of SO galaxies (see Figure [6]) which are t hought to be- 
come increasingly comm on m clusters (e.g., Dressier et al. 
|1997[ [Smith et aLl|2005 l although recent work disputes 
this ( [Holden et al.||2009 |. While we have specifically at- 
tempted to exclude SOs from our disk sample (they are 
classified as spheroidals) , distin guishing them at dis tant 
epochs remains difficult (e.g., |Moran et aL|2007b[ ). A 
more promising link exists with "passive spirals" (de- 
fined via far UV and spectra line diagnostics) which 
likely transform into SOs in cluster environments ( |Moran 



et~aLl[2005 



thought to 



2006 



2007a I . This transformation process is 
aegin slowly (1-2 Gyr timescales) in infalling 
groups, perhaps as a result of "gentle" interactions with 
other galaxies and the intergroup medium. If the intra- 
cluster medium (ICM) is sufficiently dense, ram pressure 
stripping strongly accelerates the quenching of star for- 
mation, producing SOs and ellipticals closer to the cluster 
core and dri ving evolution in the morphology-density re- 
lation (e.g., Capak et al.|2007a ). 
Because there are only 1-2 c l usters in the 1.6 deg 2 



COSMOS field flMeneux et al. 2009), cluster mecha- 



nisms cannot fully account for the passive disk popu 
lation studied here. Still, gro up scale environment s may 
play an important role (see Wilman et al. 2009). To 
investigate further, we define a "field" sample by remov- 
ing sources along the line-of-sight to dense structures. 
These were identified in t he X-ray-selected COSM OS 
group catalog ass embled by |Finoguenov et al. (20071 and 
|Leauthaud et al.| ( |2010[ ) which contains 206 groups with 
weak-lensing calibrated virial (total) masses complete to 
M 20 o ~ 3 x 10 13 M Q to z = 0.5 and M 20 o ~ 3 x 1O 14 M 
to z = 1.0. Any galaxy within the projected -R200 virial 
radius of any group, regardless of its redshift, is removed 
from the field sample, reducing the full sample by «30%. 
Remarkably, the mass dependence and evolution of pas- 
sive galaxy morphologies (Figure [3]) in this field sample 
is virtually identical to the full sample. This is not to 
say that there are no environmental trends in the pas- 
sive disk fraction which would be expected from previous 
studies of morphology-density relation. Our "field" test 
simply shows that a large fraction of passive disks live 
and evolve in halos less massive than 10 13 Mq (given our 
two lowest z-bins) , ruling out environmental processes in 
more massive halos as the sole explanation. 

Even at M200 ~ 1O 13 M environmental effects can be 
important as very weak ram pressure (2-3 orders of mag- 
nitude weaker than in cluster cores) can cause "stran- 
gulation" or "starvation" via stripping of the galaxy's 
hot gas halo, thereby shutting off a potential fuel source 
for subsequ ent star formation (see McCarthy et al.|2008 



the prevalence of passive disks at low masses, since sim- 
ulations show that lower mass galaxies are more easily 
affected. But the starved remnants should have struc- 
tural properties similar to star-forming disks, a r esult 
that seems to be at o dds with Figure [7} In addition, |Mc-| 
Carthy et al. ( |2008[ ) show that less concentrated disks 
are more susceptible to starvation. Yet, the very low- 
concentration "extreme disk" morphologies found in the 
star-forming population are almost completely absent 
among passive galaxies. Furthermore, because as M* in- 
creases starvation becomes less effective, we would expect 
the least concentrated passive disks (the "late types" ) to 
occur at the highest masses. Our observations show the 
opposite trend. We also might expect the group abun- 
dance and typical g as density, wicm. to grow with time 
from z = 1.5 (e.g.,|T^ 
massive examples of 



2009). Yet, the most 



McGee et al. 

passive disks, tor which tiicm must 
be large to induce starvation, occur in our highest red- 
shift bins (z ~ 1). 

We therefore conclude that while it is possible that 
passive disks are a product of low-mass group environ- 
ments, proposed environmental processes at these scales 
appear to have some difficulty explaining the observed 
mass dependence, evolution, and concentration of these 
galaxies. 

5.4. Gas-rich Mergers, Disk Regrowth, and AGN 
Quenching 

Beyond environmental processes, star formation 
quenching has often been associated with feedback mech 
anism s triggered by major mergers (e.g., Hopkins et al. 



2008) 
foFt 



The existence of passive disks presents a problem 
lis hypothesis: how can a major merger shut down 
star formation in these galaxies without destroying the 
underlying stellar disk? One possible solution is that be- 
fore quenching comple tes, an entirely new disk is formed 
(Hammer et al. 2005). Recent high resolutio n merger 



simulations have i llustrated exactly this effect ( Springel 



fc Hern quist 2005 ; Robe rtson et al.|2"0 06 ; Governato et a 
2009 ) . Roughly speaking, new gaseous disks can form 
m gas-rich mergers with / gas > 0.5 because a lack of 
central instabilities prevents gas transport to the center 
and suppresses bulge growth. Given observations of / gas 
as a function of M*, mergers with ps 1O 1O M are 
much more likely to result in a remnant with a disk-like 
morphology than sim ilar mergers at M» « 10 11 Mq (see 
Hopkins et aL]|2009a| ). 

One of the key aspects of the merger/regrowth model is 
that the strength of the remnant disk component dep ends 



Bekki |2009[ ) . The starvation scenario would help explain 



on / gas ( jStewart et al.|2009b| [Hopkins et aL|2009b[ ). To 
investigate how predictions trom this model compare to 
our observations, we used the relationships between / gas , 
merger mass rati o (/z), M», and remna nt bulge-to-total 
ratio {B/T) that Hopkins et al.| p009a[ ) derive based on 
a large suite of merger simulations. We assume a scatter 
of 0.25 dex in the predicted B/T and consider mergers 
with 0.3 < \i < 1.0. We use / gas (M«) evaluated at z = 2 
to correspond to passive disk progenitors at z ~ 1, and 
assume that the z = relation applies adequately to 
progenitors at z ~ 0.3. Remnant disks are those with 
B/T < 0.6. The results of this exercise are plotted in 
Figure[8]as two shaded regions for the two redshifts evalu- 
ated. These regions represent the predicted disk fractions 
spanned by mergers of different mass ratios. The region 
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Fig. 8. — Comparison of the predicted fraction of simulated 
merger remnants that would be classified as disks to the observed 
fraction of red-sequence galaxies with disk morphologies (early + 
late types). The shaded regions denote predictions for two red- 
shift ranges, z ~ 1 and z ~ 0.3. In either case, low mass ratio 
mergers predict a high abundance of remnant disks while one-to- 
one mergers trace the decreasing amount / gas as a function of M 4 . 
Overplotted are the observed passive disk fractions in the two z- 
bins. The toy model matches the observed passive disk fractions, 
as well as their mass and redshift dependence, if most formed in 
mergers with fi ~ 1.0. 



corresponding to the z = f gas {M^) is offset to lower 
M*, reflecting the decline of / gas at higher masses. Over- 
plotted are the observed fractions of passive disks (early 
+ late disks) at the two redshifts considered. While in 
detail, the scatter in predicted B/T likely depends on 
mass ratio and / gas , this toy remnant model does a good 
job matching the observed passive disk fractions as well 
as their mass and redshift dependence if most formed in 
mergers with /i ~ 1.0. A more detailed comparison be- 
tween predicted and observed morphologies is needed to 
test this further. 

This merger/regrowth scenario is consistent with the 
lack of a strong environmental dependence (Section 5.3 1 
and would predict that some passive disks will be iso- 
lated. This could be tested through clustering and pair 
count studies. Furthermore, the somewhat higher con- 
centration in passive disks compared to star-formers is 
naturally explained by the extra bulge mass added dur- 
ing the merger. The big drawback remains understand- 
ing how mergers lead to the quenching of star formation. 
If quenching is short-lived a nd violent, as in models em 
ploying a q uasar phase (e.g., Menci et al.|2006 



et al. 



Hopkins 



enough 



2008[ ), it must somehow be delayed unti. 
time passes for a stellar disk to regrow (at least 1 Gyr) 
At the same time, recent simulations have cast doubt 
on the ability of quasar-mo de AGN feedback to directly 
quench star formation (e.g., Debuhr et al.|[20To] ) . 

Instead AGN may be responsible for subtle feedback 
cycles that prevent gas cooling on longer ti mescales (1-2 
Gyr) after the cold gas is depleted (e.g., ICroton et al. 



20061. If this feedback is tied to the central black hole 



mass, this coul d ex plain why passive disks have large 
bulges (Section 5.2), but it does not make clear why 



other star-forming disk galaxies with similar bulges have 
not been quenched (Figure u\ . Finally, if quenching is 
initiated by low lu minosity AGN feedback, are mergers 



even required (e.g., Ciotti & Ostrikcr 2007)? While the 



specific quenching process in these scenarios remains un 
known, the agreement in Figure[8]ncverthcless supports a 
formation mechanism in which the frequency and appear- 
ance of passive disks is strongly linked to their progenitor 
gas fractions. 

5.5. A Two Stage Scenario for Morphological Evolution 

If star formation in passive disks is shut down be- 
cause the gas supply is depleted, simulations show that 
such gas-poor disks can be more easily disrupted by 
mergers with lower mass ratios — that is, more minor 
interactions — than disks with high / gas , owing to the 
ability o f gas to regrow disks and suppress bulge for- 
mation (Hopkins et al. 2009a). This suggests a two- 
step process where first, late-type star- forming disks are 
quenched and transformed into bulge-dominated, passive 
disks. This makes them vulnerable to minor interactions 
that, second, disrupt the disk component and convert 
them into spheroidalsp^j 

Figure[2]suggests that this two-step evolution may have 
a broad impact on interpretations of how morphological 
populations evolve. The figure shows that passive disks 
account for as much as one-third of all disk galaxies at 
M* ~ 10 11 Mq while the number declines toward lower 
masses (and redshifts). Assuming they are gas-poor and 
susceptible to more minor interactions, then because mi- 
nor mergers are far more common than major mergers 
(e.g., |Stewart et aL |2009al, this passive one-third of the 



disk population is likely to account fo r many of the newly 



formed spheroidals since z ~ 1 (e.g., Bundy et al. 2005 



Ilbert et al. 20101. But recent efforts to compare the 
formation rate of spheroidals to the merger rate have fo- 
cused only on rarer, major mergers, assuming perhaps in- 
correctly that single mergers of gas-ric h disks are the only 
channel for spheroidal formation (e.g ., Bundy et al.|2007 
Genel et al.|2008|| Bundy et al.|2009| ). These studies have 



suggested that the major merger rate, both as observed 
and predicted, may be 2-3 times too low to explain the 
formation of new spheroidals. While a more precise com- 
parison requires consistency in defining m ass ratios, mor- 
pholo gies, and understanding gas effects ( [Hopkins et al. 
2010 1 including more prevalent minor mergers involving 



passive disks could help bring current measures of the 
two phenomena into better agreement. 

5.6. Morphological Quenching 

The environmental effects and merger/regrowth model 
we have considered depend on phenomena that are ex- 
ternal to the galaxy. It is also possible that passive disks 
form as a result of internal processes . An example of thi s 
is "morphological quenching" (MQ Martig et al. 2009[ ). 
In this model, after the cold streams that fuel high- z 



star f ormation shut down or become clumpy ( Dekel et al. 
2009), the further fragmentation of gaseous disks — arid" 
therefore its ability to form stars — can be internally sup- 

18 We avoid the term "dry mergers" here, even though the inter- 
actions considered would have low / ga s, because dry mergers have 
come to denote the more narrow concept of nearly binary mergers 
between gas-poor, but already spheroidal galaxies. 
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pressed by two processes. The first is the declining self- 
gravity of the gas as it is converted to stars and the 
second is the deepening central potential caused by stel- 
lar bulge growth resulting from instabilities in the stellar 
disk or even minor mergers. The predicted remnants of 
morphological quenching are consistent with many of our 
observations. They would appear as passive disks with 
large bulges, would have a weak environmental depen- 
dence, and exhibit slow quenching, in this case initiated 
by gradual changes in internal structure. 
While more detailed comparisons between observa- 



tions and simu lations of morphological quenching ( Mar 
tig et al.]|2009 1 are clearly needed, we note that many ot 



the examples in Figures H]-[6] show large, thick disks. In 
future work, this may help (distinguish between an early 
formation with thicke ning from disk instab ilities that 



could accelerate MQ flBournaud et al. 12009 ) and later 



disk regrowth ([Hopkins et al.||2009ap . It is also impor- 
tant to determine if MO can explain why bulge growth 
leads to quenching in some galaxies but not in others 
as well as why passive disks are more disk-dominated at 
lower M„ (Figure [2]) . The higher gas fractions at these 
masses would naively require larger, not smaller bulges. 
Finally, we note that morphological quenching requires 
testa ble deviations from the Kennicutt-Schmidt relation 
(e.g., Donovan et al. 2009) and implies that gas disks 



in passive galaxies remain intact. Subsequent mergers 
of passive disks in this scenario could therefore produce 
new star-forming galaxies instead of passive spheroidals. 

5.7. Comparison with an SDSS Sample 

Finally, while this work has focused on mas s- limited 
sampl e s at h igh redshift, a detailed study by |Masters| 



et al. 



2010[ ) of a subset of passive disks identified us- 
ing GZ in SDSS provides additional insight. A direct 
comparison is not possible because the GZ passive disks 
were selected using a single color (g — r) and were re- 
quired to be mainly face-on and to have spiral arms as 
identified by GZ users. Still, the two samples paint a 
consistent picture. Red GZ disks are more common at 



higher masses (see also van der Wei et al. 2009) and 
tend to host larger bulges than their star-forming coun- 
terparts. Their emission lines indicate little ongoing SF 
but are inconsistent with post-starburst spectra, favor- 
ing instead a gradual decline of SF over 1-2 Gyr. They 
tend to live in intermediate environments consistent with 
the <1O 13 M0 halos our observations suggest. More than 
75% host nuclear bars and 50%-80% show evidence for 
LINER emission. Assuming the GZ sample is representa- 
tive of passive disks in general, the Masters et al. (20101 
analysis supports a formation mechanism that is strongly 
tied to a bulge and a central bar, somewhat weakly tied 
to environment, and capable of quenching star-formation 
slowly. 

6. SUMMARY AND CONCLUSIONS 

We have used the COSMOS survey to identify a signif- 
icant population of non-star-forming, passive disk galax- 
ies and have studied their importance with respect to 
the growth of the red sequence and the formation of 
spheroidal galaxies. After removing dusty contaminants, 
evidence for disk components is found in nearly half of 
galaxies on the red sequence at all redshifts, primarily at 
the low-mass end. The absolute number of passive disks 



below ~1O 11 M0 appears to increase with time, mirror- 
ing the growth of the red sequence, although the number 
of passive spheroidals increases faster. Above ~10 11 M Q , 
the abundance of passive disks has declined since z ~ 1 . 
We interpret these trends to indicate that passive disks 
may be a relatively common phase of galaxies transition- 
ing onto the red sequence, accounting for as much as 
60% of newly quenched systems. The declining contri- 
bution of passive disks to the red sequence with time sug- 
gests that, once formed, they transform into spheroidals 
on moderately fast (1-3 Gyr) timescales. This is likely 
aided by the fact that unlike star-forming galaxies, pas- 
sive disks may be more susceptible to morphological evo- 
lution indu ced by minor mergers , given their depleted 



gas supply (jHopkins et al. 2009b|. Thus the formation 

^Jo 



of many red sequence ellipticals does not occur in a sin- 
gle event like a major merger, and may instead proceed 
through two or more evolutionary stages. 

Understanding how passive disks form can shed light 
on the physical processes that shape galaxy evolution 
more broadly. Do passive disks evolve slowly from star- 
forming disks that exhaust their gas supply after cooling 
is suppressed or become stabilized to fragmentation? Or, 
do they form more rapidly after some event, such as a 
merger or falling into a galaxy group? While more work 
is needed, early results from this a nalysis and a study of 



a subset of passive disks in SDSS (Masters et al. 20101 
provide tantalizing clues. First, passive disks are al- 
most exclusively bulge-dominated. Star-forming disks on 
the other hand include a large number of strongly disk- 
dominated galaxies at all epochs. A bulge component is 
either required to shut down star formation — perhaps be- 
cause a bulge signals the presence of a potential AGN — or 
is built during the formation process. The presence of a 
large bulge, however, does not guarantee that star forma- 
tion is shut down. Beyond this, passive disks with the 
same mass, redshift, and morphological type are more 
concentrated on average than star-formers. While these 
assertions should be verified with stellar mass bulge/disk 
decompositions, it appears that at least some fraction of 
passive disks is structurally distinct. 

We have used a toy model to show that mergers, which 
naturally i ncrease bulge mass but can also form a new 
disk (e.g., |Springel fc Hernquist 20051, produce rem- 
nants with a similar mass and redshift dependence as 
passive disks. This agreement suggests a strong link 
between morphology and / gas - The problem with the 
merger/regrowth model remains understanding how the 
merger induces quenching. SDSS observa tions of passive 
disk emission lines favor a slow process (Masters et al. 
20101 that operates over 1-2 Gyr. This seems inconsis- 



tent with a merger-triggered quasar phase, which must 
be delayed to allow enough time for a new disk to form. 
On the other hand, slow quenching may take place in 
1O 13 M halos in which the group medium can slowly 
"starve" galaxies of fuel by stripping their hot gas reser- 
voirs. Environmental processes therefore present an al- 
ternative to the merger/regrowth model but have trou- 
ble accounting for the importance of bulges and the way 
in which passive disk morphologies depend on mass and 
redshift. The newly propose d morphological quenching 
scenario (Martig et al. 2009) also predicts slow quench- 
ing as internal structural changes suppress disk frag- 
mentation and star formation. This mechanism requires 
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large bulges, as observe d, and predicts thick stellar disks 
More detailed comparisons and 



(Bournaud et al 



2009) 



a test oi the predicted mass and redshift dependences in 
this scenario are now needed. 

The prospects for resolving the formation mechanism 
look promising and we note that some combination of 
the processes considered here may be operating simulta- 
neously. More detailed studies of the structure, color gra- 
dients, and stellar populations of passive disks will help 
link them with their progenitors. At the same time, it is 
important to verify the high disk fractions revealed by the 
ZEST classifier with bulge-to-disk decompositions, ide- 
ally performed in terms of stellar mass. This will confirm 
the significance of disk components and enable more di- 
rect comparisons to star-forming samples as well as simu- 
lations. Improvements at low- 2: will come from expanded, 
mass-limited samples, and at high-z from clustering and 
pair-count studies. Larger samples from future surveys 
will overcome cosmic variance and allow individual pop- 
ulations to be tracked as they evolve through different 



stages. Both with current surveys and in the future, 
passive disks offer a unique opportunity to separate and 
understand the processes that drive global patterns of 
galaxy evolution since z ~ 1. 
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